Test for interlayer coherence in a quasi-two-dimensional superconductor. 
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Peaks in the magnetoresistivity of the layered superconductor «:-(BEDT-TTF)2Cu(NCS)2, mea- 
sured in fields < 45 T applied within the layers, show that the Fermi surface is extended in the 
interlayer direction and enable the interlayer transfer integral {t± ~ 0.04 meV) to be deduced. 
However, the quasiparticle scattering rate is such that h/r ^ 6t±, implying that k-(BEDT- 
TTF)2Cu(NCS)2 meets the criterion used to identify interlayer incoherence. The applicability of 
this criterion to anisotropic materials is thus shown to be questionable. 
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Many correlated-electron systems which are of funda- 
mental interest have very anisotropic electronic band- 
structure. Examples include the "high-Tc" cuprates 
H, layered ruthenates and crystalline organic met- 
als 1^, |j. Such systems may be described by a tight- 
binding Hamiltonian in which the ratio of the interlayer 
transfer integral t±^ to the average intralayer transfer inte- 
gral i|| is < 1 ||, |, H. The inequality H/t >tj_^ where 
is the quasiparticle scattering rate ^, ||] , frequently 
applies to such systems, suggesting that the quasiparti- 
cles scatter more frequently than they tunnel between 
layers. The question has thus arisen as to whether the 
interlayer charge transfer is coherent or incoherent, i.e. 
whether or not the Fermi surface (FS) extends in the in- 
terlayer direction H, ||, ||. In this paper we have used 
magnetoresistance data to estimate the interlayer trans- 
fer integral in the highly anisotropic organic supercon- 
ductor K-(BEDT-TTF)2Cu(NCS)2. We find that the ma- 
terial obeys the inequality h/r > tj_; moreover, mean- free 
path in the interlayer direction is <20% of the unit-cell 
height. Nevertheless, our data demonstrate a FS which 
is extended in the interlayer direction. 

K-(BEDT-TTF)2Cu(NCS)2 was selected for our exper- 
iments because it is perhaps the most thoroughly char- 
acterised quasi- two-dimensional (Q2D) conductor Q. In 
contrast to the cuprates, the FS topology is well known 
from Shubnikov-de Haas (SdH) and de Haas-van Alphen 
(dHvA) studies Q and from angle-dependent magne- 
toresistance oscillation (AMRO) |^ and millimetre-wave 
(MMW) experiments it consists of a pair of quasi- 
one-dimensional (QID) electron sheets plus a Q2D hole 
pocket (see Fig. |f|a [|lO|, 0). The K-phase BEDT-TTF 
salts are considered to be leading contenders for inter- 
layer incoherence and optical data may be inter- 
preted as consistent with this suggestion ||l2[. More- 
over, models for unconventional superconductivity in bi- 
phase BEDT-TTF salts invoke the nesting properties 
of the FS 0, |I3|, the degree of nesting might de- 
pend on whether the FS is a 2D or 3D entity (see [Q, 
Section 3.5). Experimental tests for coherence in n- 
(BEDT-TTF)2Cu(NCS)2 are thus far deemed to be in- 



conclusive g; e.g. semiclassical models can reproduce 
AMRO § and MMW data § equally well when the in- 
terlayer transport is coherent or "weakly coherent" [^. 

To examine how interlayer coherence might be de- 
tected, we use a tight-binding dispersion relationship 
~2t± cos(fc^a) to represent the interlayer dispersion, 
where kj_ is the interlayer component of k and a is the 
interlayer spacing. This is added to the effective dimer 
model, which is known to represent the intralayer band- 
structure accurately lIC 
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to yield 



£;(k) = ±2cos{^)^tl^+tl2 + 2t^it^2 cos(fccc) 



-|-2tb cos(fcbfe) — 2t± cos{k±a). 



(1) 



Here kh and fee are the intralayer components of k (see 
Fig. nla) and th, tci and tc2 are interdimer transfer inte- 
grals ||ig, |ll| ; the + and — signs result in the QID sheets 
and Q2D pocket of the FS respectively (Fig. |l|a). The 
addition of the interlayer dispersion produces a warp- 
ing of the FS, shown schematically for the Q2D section 
m Fig. |l|b; the FS cross-section is modulated in the in- 
terlayer direction. This modulation might suggest that 
separate "neck and belly" frequencies would be observed 
in the dHvA effect however, only a single frequency 
is seen in low-field studies iQ, suggesting that the cy- 
clotron energy exceeds t± at the fields at which quantum 
oscillations are observed. 

Having failed to detect the beating between neck and 
bellies, and in view of the inconclusive nature of other 
tests we have chosen instead to examine the be- 
haviour of K-(BEDT-TTF)2Cu(NCS)2 in almost exactly 
in-plane magnetic fields |||, |T^, |f8| . The motion of quasi- 
particles of charge g in a magnetic field B is determined 
by the Lorentz force Kk — gv x B, where the quasi- 
particle velocity is given by Tiv = Vk-B ||l^, 0; this 
leads to orbits on the FS in planes perpendicular to B. 
Hence, if the FS is extended in the interlayer direction, 
an in-plane field will cause closed orbits on the bellies 
(see Fig. ^d). Such orbits are effective at averaging v±, 
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the interlayer component of the velocity, and their pres- 
ence will lead to an increase in the magnetoresistivity 
component p^z 0, B can then be tilted away 

from the in-plane direction by an angle A, such that the 
small closed orbits about the bellies cease to be possible; 
this occurs when B is parallel to v at the point at which 
is a maximum (i.e. when the normal to the FS is 
at a maximum angle to the Q2D planes). Therefore, on 
tilting B around the in-plane orientation, we expect to 
see a peak in p^z, of angular width 2 A, if (and only if ||]) 
the FS is extended in the interlayer direction. 

A problem in using k-(BEDT-TTF)2Cu(NCS)2 is its 
high in-plane critical field; PqHc2{T = 0) « 35 T, falling 
to ijlqHc2 w 25 T at 4.2 K ||21|. Moreover, apparent 
peaks in the resistivity of k-(BEDT-TTF)2Cu(NCS)2 in 
in-plane fields may occur due to dissipative processes as- 
sociated with vortices within the mixed state ^ . To 
ensure that such effects do not interfere with our data, we 
choose to stay at fields well above the superconducting 
state; our measurements were therefore carried out in the 
45 T hybrid magnet at NHMFL Tallahassee. The experi- 
ments involved two single crystals (~ 0.7 x 0.5 x 0.1 mm'^; 
mosaic spread <0.1°) of k-(BEDT-TTF)2Cu(NCS)2, pro- 
duced using electrocrystallization . In one of the crys- 
tals, the terminal hydrogens of the BEDT-TTF molecules 
were isotopically substituted by deuterium; we refer to 
the deuterated sample as d8, and the hydrogenated sam- 
ple as h8 |2^. Both crystals were mounted in a ^He cryo- 
stat which allowed rotation to all possible orientations in 
B [Q; sample orientation is defined by the angle 9 be- 
tween B and the normal to the sample's Q2D planes and 
the azimuthal angle (/)(</) = is a plane of rotation of B 
containing kb and the normal to the Q2D planes). The 
interlayer magnetoresistance Rzz (oc Pzz) was measured 
using the four-terminal ac techniques described in 

Fig. ^ shows Rzz of the d8 sample close to the in- 
plane orientation = 90°; data for three values of 4> are 
shown. (The data for the hS sample were very similar 
in all respects [^.) The edges of Fig. |^ are dominated 
by AMROs and related phenomena; as these are well 
known |, H in «;-(BEDT-TTF)2Cu(NCS)2 g, we 
shall not describe them further in this paper. Close to 
— 90°, there is a distinct peak in Rzz, the width and 
height of which vary with 0; we attribute this peak to 
the closed orbits described above. 

Fig. H shows the peak at temperatures T from 0.48 K to 
5.1 K; increasing T by over an order of magnitude reduces 
the peak height but has little effect on its width [^ . 
Varying the field (in the range 35 — 45 T at 0.5 K, and in 
the range 29 — 45 T at 4.2 K) also has little effect on the 
peak width; increases of field result in increases in peak 
height and definition, in a manner similar to the field 
dependence of AMROs |l^. Both of these observations 
support the idea that the peak is a consequence of the 
FS geometry alone. 

Fig. ^ shows the variation of 2A, the full width of 
the peak close io = 90° versus 0; the width was de- 
duced using the extrapolations shown Fig. (inset). In 



order to interpret this variation, we use Eqn. |lj to cal- 
culate w_Lmax, the maximum value of v^, and v\\, the 
intralayer velocity component parallel to the plane of ro- 
tation of B; when measured in radians, A « f^max/fn. 
As all of the relevant quasiparticle motion occurs close 
to the Fermi energy, E w iJp, we adjust the parameters 
of Eqn. to reproduce the known FS parameters of k- 
(BEDT-TTF)2Cu(NCS)2. First, (where U is the 

mean of tci and tc2) is adjusted to obtain the dHvA fre- 
quencies of the Q2D pocket and the magnetic breakdown 
orbit [^0|, |2^, The absolute value of tc is then con- 
strained by fitting to the effective mass of the breakdown 
or bit |, |2|. Third, the energy gap measured in mag- 
netic breakdown [E^ « 7.8 meV [^) gives tci — ic2 = 
Eg/2 0, leading to tb = 15.6 meV, tci = 24.2 meV and 
tc2 = 20.3 meV. Finally, a = 16.2 A ||T|], so that equa- 
tions for A contain only one adjustable parameter, t±. 
The substitution of t± = 0.04 meV leads to the curves 
shown in Fig. ^; the continuous curve is from the Q2D 
FS, and the loops are caused by the QID sheets, which 
support closed orbits only over a restricted range of (j). 
In the latter case, the lower branch of the loop results 
from the flatter portions of the QID FS, roughly parallel 
to kc, whereas the upper branch is caused by the more 
highly curved region at the zone boundary (Fig. |^a). 

For most (f), there is agreement between curve and data 
(Fig. ^), but around cf) = and 180° it seems that the 
observed peak width is sometimes dominated by closed 
orbits on the QID sheets, and sometimes by those on 
the Q2D FS section; the dominant width presumably de- 
pends on which FS section is more effective at averaging 
v± ||l^ . MMW studies |^ suggest that the interlayer cor- 
rugations of the QID sheets are more complex than those 
given by Eqn. |l|. This could lead to a rapid variation with 
(j> of the effectiveness of the orbits in increasing pzz ■ For 
a narrow range of angles close to (j> — 25°, the peak was 
particularly wide and large (see top trace. Fig. This 
is perhaps connected to the gap between the QID and 
Q2D FS sections. 

As mentioned above, the bandstructure of k-(BEDT- 
TTF)2Cu(NCS)2 within the Q2D layers is determined 
by the interdimer transfer integrals |ll|. A better 
guide to the total intralayer bandwidth than the param- 
eters used above (which are relevant for E w Ep and 
thus include renormalising interactions 0) is given 
by the fits to optical data of Ref. [0, which suggest 
tb ~ 60 meV and tc « 120 meV; these are a factor >10'^ 
larger than t± « 0.04 meV. 

The failure of the dHvA effect |l^ to observe necks and 
bellies may now be understood; the Landau-level spacing 
at the lowest |B| used (- 6 T) is 0.2 meV, 5tj_. 

It is instructive to compare t± with r^^. Samples 
d8 and h8 have been studied using SdH oscillations at 
|B| < 15 T At such fields, the oscillatory magnetore- 
sistance is much less than the nonoscillatory component, 
and magnetic breakdown is a minor consideration [Q. 
Hence, the 2D form of the Lifshitz-Kosevich formula may 
be used to extract r giving r = 2.9 ± 0.5 ps (h8) and 
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T = 2.6 ± 0.3 ps (d8) Another estimate can be 

derived from MMW studies which measure the FS- 
traversal resonance (FTR) due to quasiparticles crossing 
the QID FS sheets; these experiments used samples from 
the same batch as h8. In the data of the FTR ap- 
pears at i? « 10 T, with a full-width at half-maximum 
of AB « 7 T. If we assume that ujt - B/AB where 
w = 27r X 70 X 10^ rad s~^ we obtain r ^ 3 ps, close 
to the SdH values. Thus, H/t « 0.24 meV, ~ 6t±. 

Finally, we estimate the intralayer mean-free path 
A||; a typical value of W|| (w 7.2 x 10'' ms~^), yields 
A|| = viiT « 0.2 /zm {i.e. k-(BEDT-TTF)2Cu(NCS)2 is 
a very clean system). However, the maximum interlayer 
velocity is wj^max ~ 2t±a/h « 200 ms^^, suggesting an 
average interlayer velocity v± ^ 100 ms^^, and a mean- 
free path A_L ~ v±T w 3 A. The usual mean-free path for- 
malism implies that only ^ exp(— 16.2/3) w 0.5% of the 
quasiparticles travels between adjacent layers (separation 
16.2 A [|l5|) without scattering. In such circumstances, 
conventional criteria Q lead one to expect incoherent in- 
terlayer transport, yet the peak in p^z described above 
unambiguously demonstrates a 3D FS topology. 

The criterion for incoherent transport was developed 
for isotropic disordered metals [Q; the mean free path A 
was taken to represent the spatial extent over which the 
Bloch waves are coherent and is, by default, isotropic. 
However, in Q2D systems A|| ^ A_l, so that a purely 
ballistic model of transport (of the sort used to treat dis- 
orded 3D systems Q]) fails to preserve information on 
the spatial extent of the wavefunctions perpendicular to 



the layers. Our experiment therefore demonstrates seri- 
ous inadequacies in such models. When impurities and 
defects are randomly dispersed throughout a material, 
the inter-impurity separation should be almost isotropic. 
Hence, we should also expect the spatial coherence of the 
Bloch states to remain roughly isotropic, even when the 
group velocity is highly anisotropic. 

In summary, we observe a peak in the interlayer resis- 
tance of the highly anisotropic superconductor k-(BEDT- 
TTF)2Cu(NCS)2 when a magnetic field is applied within 
the layers. This demonstrates Q that the Fermi surface 
is extended in the interlayer direction, and allows the 
interlayer transfer integral to be estimated to be t± ~ 
0.04 meV {t±/kB ~ 0.5 K). The small size of t± explains 
unsuccessful attempts to observe neck and belly orbits 
in dHvA studies. However, perhaps the most interesting 
point to emerge is that k-(BEDT-TTF)2Cu(NCS)2 obeys 
the criteria commonly used to delineate interlayer inco- 
herence - the interlayer mean free path is smaller than 
the unit cell size, and h/r > t± - and yet it clearly pos- 
sesses a three-dimensional Fermi surface. It is perhaps 
now time to re-examine the validity of such criteria when 
applied to strongly anisotropic systems. 
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FIG. 1: (a) Cross-section of the Fermi surface (FS) and 
Brillouin zone of k-(BEDT-TTF)2Cu(NCS)2 predicted by 
Eqn. |l| 0. (b) Perspective view of the Q2D FS section 
described by Eqn hi; the intralayer curvature and interplane 
warping have been exaggerated for clarity. The lines indicate 
quasiparticle orbits on the FS due to the in-plane field B. 
Note the closed orbits about the "belly" of the FS. 
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FIG. 2: Interlayer resistance R^z for the d8 k;-(BEDT- 
TTF)2Cu(NCS)2 sample as a function of tilt angle 6. Data 
for three planes of rotation of the field are shown, (j> = 25° 

(upper), (j) ~ 20° (middle) and (f> = 15° (lower). The static 
magnetic field is 42 T, and the temperature is 520 mK. 
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FIG. 3: Interlayer resistance (Rzz) versus angle 9 for temper- 
atures r = 0.48 K, 1.4 K, 3.0 K, 4.4 K and 5.1 K (</> = 135°). 
The background magnetoresistance increases with increasing 
T, whereas the peak at ^ = 90° becomes smaller. The data 
shown are for the d8 k-(BEDT-TTF)2Cu(NCS)2 sample. The 
inset shows the intersections of the linear extrapolations used 
to determine the peak width. 




FIG. 4: Angular width 2A of peak in iZ^^ (B = 42 T, T « 
500 mK) versus azimuthal angle (f>; data for the d8 k-(BEDT- 
TTF)2Cu(NCS)2 sample are shown. Points are data; the 

curves represent the model prediction with t± = 0.04 mcV. 
The continuous curve is due to the Q2D FS section; the QID 
sheets can only support closed orbits over a restricted range 
of 0, leading to the top-shaped loops. 



